Hyperfine-induced decoherence in triangular spin-cluster qubits 
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We investigate hyperfine-induced decoherence in a triangular spin-cluster for different qubit en- 
codings. Electrically controllable eigenstates of spin chirality (C z ) show decoherence times that 
approach milliseconds, two orders of magnitude longer than those estimated for the eigenstates of 
the total spin projection (S z ) and of the partial spin sum (S12). The robustness of chirality is due 
to its decoupling from both the total- and individual-spin components in the cluster. This results 
in a suppression of the effective interaction between C z and the nuclear spin bath. 
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Introduction - Molecular nanomagnets represent a 
varied class of spin clusters, whose physical properties 
can be extensively engineered by chemical synthesis [![. 
This makes them candidate systems for the implemen- 
tation of spin-cluster qubits While most of the 
attention has been so far focused on the use of the total- 
spin projection (S z ) as a computational degree of freedom 
(DOF), it has been recently realized that alternative en- 
codings would enable the use of electric - rather than 
magnetic - fields for the qubit manipulation [f|. In par- 
ticular, transitions between states of opposite spin chi- 
rality [C z = (4/v3)si • S2 x S3] can be induced in antifcr- 
romagnetic triangles with Dzyaloshinskii-Moriya interac- 
tion. Spin-electric coupling constants compatible with 
ns gating times T g have been predicted by effective mod- 
els [5|, @ and microscopic ab initio calculations Q . Fur- 
ther investigation is indeed required in order identify spe- 
cific molecular nanomagnets with large spin-electric cou- 
pling, or to enhance such coupling by introducing suitable 
chemical substitutions in existing systems 0. 

In order to assess the suitability of spin chirality for 
applications in quantum-information processing, its r g 
has to be contrasted with a characteristic decoherence 
time Td- At low temperatures, quantum coherence in 
molecular nanomagnets is limited by the coupling to the 
nuclear spin environment, with typical values of Td in the 
microsecond range [9l-ITTj . All the existing literature is 
however concerned with linear superpositions of different 
S z eigenstates. Here we theoretically investigate the de- 
pendence of hyperfine-induced decoherence on the qubit 
encoding within a prototypical spin-cluster qubit, con- 
sisting of an antiferromagnetic spin triangle. In particu- 
lar, we consider three different DOF, namely S z , C z , and 
the partial spin sum S12 (S12 = Si + S2), whose value 
- like that of C z - can be controlled through the spin- 
electric coupling. Since the optimal candidate system 
has not been identified yet, we refer here to a prototypical 
molecular spin-cluster qubit, with a typical electron-spin 
Hamiltonian [l2[ and bath of nuclear spins (l3j . While 
the quantities of interest might to some extent vary from 
one molecular nanomagnet to another, the hyperfine- 
induced decoherence presents striking differences in the 
three considered DOF, that are not expected to depend 
on the specific features of the spin-cluster qubit. 

Qubit encodings in the spin triangle — We consider a 
triangle of s = 1/2 spins, with dominant antiferromag- 



FIG. 1: (color online) (a) Schematics of the local spin pro- 
jections (sz t i) in the spin triangle, corresponding to the |0) 
(red) and |1) (blue) states in the three considered qubit en- 
codings. All the states have S — 1/2; besides, {k\S z \k)c z = 
(k\S z \k) Sl2 = -1/2 and {k\C z \k) Sz = 1, with k = 0, 1. "(b) 
Angle ak between the vector (S) and z for the eigenstates 
|0)c, (blue) and \l)% z (green), for A/gfi B = 0.5 (solid lines) 
and 2.0 (dotted), (c) Statistical distribution of the distances 
d e n between the A^ e = 3 electron and the N„ — 200 nuclear 
spins with randomly generated positions. 



netic coupling and Zeeman interaction: 



H = J^ Si 



3/i S B • S. 



(1) 



An additional term Hi determines the expression of the 
lowest eigenstates |0) and |1), belonging to the ground 
state 5 = 1/2 quadruplet. As discussed in the following, 
the robustness of the spin-cluster qubit with respect to 
hyperfine-induced decoherence strongly depends on the 
distinguishability between |0) and |1) in terms not only 
of total spin orientation, but also of spin texture. Here- 
after, we thus discuss these features in some detail in two 
relevant cases: 



Bf- = dz-y: 



H? 12 = (J12- J) 81 -82. 



(2) 
(3) 
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FIG. 2: (color online) Time dependence of the decoherence 
factor r for the three qubit encodings: S z (black), 512 (red), 
and C z (green for 8 — and blue for 9 = 7r/8). The curves 
are averaged over Nx — 5 x 10 4 randomly generated initial 
states \T) of the nuclear bath. Inset: Statistical distribution 
(squares) of the parameter Si, and corresponding Gaussian 
fits (solid lines); same convention as above for the colors. 



The term H 1 * accounts for the Dzyaloshinskii-Moriya 
interaction in a spin triangle with, e.g., D%h symmetry 
12j . For H e = H + , the four lowest eigenstates 
can be labelled after the value of the spin chirality C z — 
(4/ -\/3)si -S2 x S3, and the Dzyaloshinskii-Moriya term can 
be rephrased as: H?* = AC Z S Z , with A = Dy/?> @. In 
particular, if the magnetic field is oriented parallel to the 
principal axis of the molecule (B || z), the eigenstates 
\C Z ,S Z ) read: | ± 1, +1/2) - (| Iff) + e ±l2 ^ 3 | tit) + 
eT i2./3| m) )/V3 and I ± 1, -1/2) = a x a x a x \ ± 1, +1/2), 



where a x is the Pauli operator acting on Si. Both S z 
and C z commute with the electron-spin Hamiltonian H e , 
which makes them suitable as computational DOF. In 
the first case, the logical states are: 

|0>s = |S^-1/2;C Z =+1), |l) s = |S»=fl/2; C z =+l). 

The expectation values of are oriented along the mag- 
netic field, are identical for the three spins, and change 
sign with the qubit state [Fig. [T](a)]: 



(lK<|l}s. = -<0|«,,i|0) s . = 1/6. 



(4) 



If the computational DOF is identified with spin chirality, 
the logical states are instead: 

|0)c,= |C^+l;^=-l/2), \l) c = \C^-l;S z =-l/2), 

and the expectation values of the three spins are inde- 
pendent on the qubit state [Fig. Q] (a)]: 



(1K<|1>C = (0|a»,i|0) o , = -1/6. 



(5) 



Such condition is however not general. In fact, if the 
applied magnetic field is tilted with respect to the z axis, 
B = £>(sin6>x + cos6>z), C z is still a good quantum num- 
ber, but (k\si\k)c (with k = 0, 1) are always oriented 



along B J. = (B x ,0,B z ±A/g/j,B)- Eigenstates of opposite 
chirality are thus characterized by different orientations 
of the spin expectation values [see Fig. QJb)]: 

(M 8 x,i\k)c z = sinafe/6, (k\s Zii \k) e Cz = cosa fc /6, (6) 



where ctk = arctan 



\B sin ( 



-7T, < arctan < 



B cos8+(-l) k A/gfj, E 

7r and x = ±1 for A ^ B gus- 
li no Dzyaloshinskii-Moriya interaction is present and 
one exchange coupling differs from the other two, the 
term iff* is replaced by iff 12 (Eq. EJ. For H e 
If, + iff 12 , the four lowest eigenstates can be labelled 
after the partial sum of the first two spins, rather than 
the spin chirality: \Si2,S z ), where S12 — 0,1. Their 
expressions read: |0, +1/2) = (| |ff) - I lff))/V2, 

ii,+i/2) = (i nt> + 1 iff) - 21 m))/Ve, while 

[512,-1/2) = alal<7 x \S 12 ,+l/2). Choosing S 12 as the 
computational DOF, one has: 



|0) Sl2 = \Sia = 0;5 2 =-l/2), |l) Sl = \S 12 = l;S z =-l/2). 

As far as the spin expectation values are concerned, S12 
represents an intermediate case between S z and C z . The 
qubit states have in fact identical values for the total 
spin, (0|S|0)s 12 = (l|S|l)s 12 , like C z , but they strongly 
differ in terms of spin texture, like S z [Fig. [Ha)] 



(0Ki=i, 2 |0) Sl2 =0, (C 
{l|**,i=l,2|l)fiT„ = -1/3, 



K 3 |0)s 12 = -1/2 (7a) 
<iK 3 |i) Sl2 = 1/6. (7b) 



Nuclear spin and hyperfine interactions — The deco- 
herence of the spin-cluster qubit is investigated by simu- 
lating the coupled dynamics of electron and nuclear spins, 
induced by the Hamiltonian H = H e + H n + H en . The 
qubit and the nuclear environment are initialized respec- 
tively in the linear superposition |-0 e (0)) = ^= (|0) + |1)) 

and in the mixed state p n {0) — ^ I Pi|I)(X|. Here, the 
expressions of |0) and |1) depend on Hi (see above), 



while |Z) 



|mf,. 



. 7)1 



N n / 



and mf are the projec- 



tions along the magnetic field direction of the N n nu- 
clear spins. In the pure-dephasing regime, each state 
|*i(0)) = ^3 (|0) + |1)) ® \X), evolves into: \* x (t)) = 

■j=(\0,Xo) + where \2 ) can be regarded 

as the state of the nuclear bath conditioned upon the 
qubit being in the |0) (|1)) state. The degree of co- 
herence in the reduced density matrix of the qubit, 
Pe = Tr„{^ I Pi|*i(t))(*i(i)|}, is given by the so- 
called decoherence factor: r(t) = Y2i Px r x(t), with 
r T {t) = (Zi(t)|2b(t)) and <0|p e |l> = rx/2. 

The nuclear spin bath we consider consists of N n = 200 
hydrogens (I = 1/2), whose positions r™ are randomly 
generated so as to reproduce typical values of the spin 
density and the electron- nuclear distances d en = |rf — r™|, 
where rf are the positions of electron spins [Fig. [1] (c)] 
(l4T |. The nuclear-spin Hamiltonian H n includes Zee- 
man and dipole-dipole terms: H n = B • ^ p lu p 1 p + 

D nn E p <J I p ' !« - 3 ( I p ' *W)(I<z • fpg)]Apg, where D nn = 
(/io/47r)/x^7| and r pq = r'' p l — r£. Electron and nu- 
clear spins are coupled by dipole-dipole and contact in- 
teractions: H en = D en Y,iJ2 p [si ■ I P - 3(sj • T ip )(Ip ■ 
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*ip)]/ r ip + J2i a i S i ' Ig(i)i wnere D en = (MO / ' ^) A* 

and r^p = r^ — r™. The contact terms a,, whose effect will 
be considered in the final part of the paper, couples elec- 
tron and nuclear spins belonging to the same magnetic 
center. 

The dephasing arises from the qubit-state dependent 
dynamics of the nuclear bath, generated by an effective 
Hamiltonian T-L. We derive T-L from the above specified 
H = H e + H n + H en in a two-step procedure [H, . 
We first project the single-electron spin operators s a ^ 
onto the 5 = 1/2 quadruplet: Vs=i/2 s p'Ps=l/2 = 

J2i,j=o( i \ s p\j) a 'ij' wnere Vij = K)0'l and K) are tne 
eigenstates of H e . We then apply a Schrieffer- Wolff 

transformation, that removes from the Hamiltonian the 
terms that are off- diag onal in the basis of electron-spin 
eigenstates \i) [l7], Il8j , and finally neglect energy non- 
conserving terms (secular approximation). The resulting 
Hamiltonian reads: H — J2k=o 1 ® Hk, where 

p=l p^q 

and z' = H/B. Differences between Hq and Hi 
result from the hyperfine interactions, and are re- 
sponsible for the qubit decoherence, being rx(t) — 
(I\exp{i'Hit/h} exp{—i'Hot/fi}\Z)- In particular, the 
quantities (wp 1 — ) are linear in H en , and essentially re- 
sult from differences in the magnetic field induced by the 
nuclear spin I p at the electron-spin positions r\ (see be- 
low). The terms (A° g — A pg ) and (B® — B pq ) are instead 
quadratic in H en , and result from the qubit-state depen- 
dence of the couplings between pairs of nuclei, mediated 
by virtual transitions of the electron spins. The time 
evolution of the nuclear states \Xk) is computed within 
the pair-correlation approximation, where the nuclear dy- 
namics is traced back to independent flip-flop transitions 
between pairs of nuclear spins jl7l [l9j . 

Hyperfine-induced decoherence - The fastest con- 
tribution to dephasing in the spin-cluster qubit is re- 
lated to inhomogeneous broadening, and typically takes 
place on time scales that are much shorter than those 
characterizing the dynamics of the nuclear bath (r„ ~ 
ft/jB* g | ~ 10 2 /is). Such contribution results from the 
renormalization of the energy gap between the states 
|0) and |1) induced by the hyperfine interactions: 8% — 

Efe=o,i(-l) fc (^^l^9.- E ? «-<K- Being 
the nuclear spin bath initially in a mixture of states 
|2), the decoherence factor evolves as: r(t -C r„) ~ 
e - i (B _B 1 )i^ 2 .p :re -i5 I t ) while \x k (t < T „)) ~ \X). In 

first order in H en , Sx can be regarded as a function of 
the Overhauser field at the electron-spin sites: 

Sx c ^g^lM) • [(0N0) - (118*11)], (9) 

i 

where Bf^rf) = £> e „ £ p m^z' - 3(z' • r tp )i lp }/rf p . 
In the case of the S z qubit (see Eq. H]), S x z ~ 
— (l-t-Bg/3)J2i-Bhfz>( r i)- The statistical distribution 
N(S X z ) is reported in the inset of Fig. [5] (black squares) 
for 5 x 10 4 initial nuclear states \X), randomly gener- 
ated from a flat probability distribution Px = l/2 Nn . 
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FIG. 3: (color online) Time evolution of the decoherence fac- 
tor r m in the cases of the Sz (black) and S12 (red) DOF. The 
case of chirality is displayed for — (green) and for finite 
values of the tilting angle (blue). All curves are averaged 
over 10 2 randomly generated initial conditions \T); the spin 
Hamiltonian parameters are: A = 1 K, B — 1 T. 



N(6 X ') is well fitted by a Gaussian function (solid line), 
with as z = 9.0 neV. Correspondingly, the decay of \r(t)\ 
(black line in Fig. [2]) is approximately Gaussian, and its 
characteristic time scale is 10 2 ns. In the case of the S12 
qubit, the three electron spins are no longer equivalent: 
4" ~ -(/W3)[2B^,(r|) - Bl f z ,{r\) - S^ z ,(r|)]. 

However, the statistical distribution of <5j 12 strongly re- 
sembles that of S z (see the red squares in the figure inset, 
and the Gaussian fit with as 12 — 9.4 neV), and so does 
the time evolution of the decoherence factor (red curve 
in the main panel). In fact, since the distances d ee be- 
tween electron spins are larger than the smallest d en [see 
Fig- Ec)] 03 > the spatial fluctuations of the Overhauser 
field within the spin triangle are comparable to its aver- 
age value. In spin clusters with larger d en /d ee ratios (not 
shown here), spatial fluctuations of B^/(r) are relatively 
small. As a result, <5f 12 <C <5f z , and the S12 qubit is less 
affected by inhomogeneous broadening than S z . 

In the case of the C z qubit and for B || z, the Over- 
hauser field does not renormalize the energy difference 
between the states |0) and |1), that have identical ex- 
pectation values for all single-spin projections (Eqs. [5J 
[5]). The leading contribution to is given by terms 
that are second order in the hyperfine Hamiltonian, 
S x z = J2 P ^q(Apq ~ Ap^rn^m^ , and its modulus is here 
5 orders of magnitude smaller than that of S x z and S x 12 . 
Correspondingly, no inhomogeneous broadening occurs 
in the considered time scale (green curve). For a tilted 
magnetic field (9 ^ 0), states of opposite chirality have 
different expectation values (sj) (see Eq. [5]), and thus 
couple differently to the Overhauser field. The leading 
contribution to the renormalization of the energy differ- 
ence reads: 8% z (6) ~ (VBg/6)J2l=o(- 1 ) k ( s ' m akB% f x ,+ 
cos akB^f z ,), where x' _L z' and lies in the xz plane. The 
statistical distribution of 5 X Z {9 — t/8) and the resulting 
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FIG. 4: (color online) Time dependence of the decoherence 
factor in the presence of three additional nuclear spins (N n — 
203) localized at the electron-spin positions rf and coupled 
to the respective electron spins with a contact coupling a p = 
lOmK (panel a), or a p = lmK (panel b). The solid curves 
correspond to the cases S z (black), S12 (red) and C z (figure 
insets) . The dotted lines represent the time dependence of r m 
in the absence of the three nuclei with contact couplings and 
correspond to those displayed in Fig. [2] 

qubit dephasing are reported in Fig. [5] (ac z — 4.5 neV, 
blue curve). 

Nuclear spin evolution in the field of electron spins cre- 
ates entanglement in the initially separable state |5'i(0)), 
and causes additional decoherence on the time scale of 
nuclear dynamics, r„. In order to single out this contri- 
bution, we compute the function: r m {t) = J2x Pi\ r i(t)\- 
In the case of the S z qubit, electron-nuclear correlations 
result in a decay of r m in the /is time scale (Fig. |H 
black curve). The decay is induced by the interplay of 
the dipolar interactions between the nuclei and of the 
term J^pi^p ~ whose expectation value gives 5 X Z 

[20| . A similar time dependence for r rn is obtained in the 
case of the S12 qubit (red curve). Here, the same terms 
in the effective Hamiltonian % dominate, and have sim- 
ilar expectation values: S x 2 ~ S x 12 (see the inset of Fig. 
[2]). This quantity (5 X *) is about 5 orders of magnitude 
smaller for the C z qubit, if B || z. As a result, the dy- 
namics of the nuclear bath is largely independent on the 
qubit state in the considered time range, and no appre- 
ciable decoherence takes place (green curve). This is no 
longer the case if the magnetic field is not aligned with 
the principal axis of the spin triangle: in fact, the deco- 
herence time rapidly decreases as S x z (8) increases with 
the tilting angle 9 (blue curves). 

We finally investigate the possible contribution to de- 
coherence of the contact terms. Such contribution is 
qualitatively different from that of the dipolar interac- 



tions, for it results from the relatively strong coupling 
with the electron spins of few (iVJj ~ N e <C N n ) nuclei. 
Here, the N% = N e = 3 additional nuclear spins are local- 
ized at the electron spin sites r|, and are assumed for sim- 
plicity identical to the remaining N% — N n — N% = 200 
nuclei. The inequivalence between the N% and N% nu- 
clear spins, resulting from strong coupling of the former 
ones with the electron spins, warrants the factorization of 
the decoherence factor: r(t) = r c (t)rd(t) [21]. The time 
evolution of r c (t) is reported in Fig. 01 for a p = 1 mK 
and a p = 10 mK [panels (a) and (b), respectively]. In the 
case of S z (black curve) and Si 2 (red), r c is responsible 
for the fast oscillations, while the decay is due to (dot- 
ted lines). Oscillations of the decoherence factor caused 
by the contact interactions are also present in the case 
of C z qubit (figure insets). These take place on a time 
scale which is much longer than that of S 2 and S12, but 
much shorter than the one that characterizes the contri- 
bution to the decoherence of C z due to dipolar coupling. 
The chirality qubit also presents a different dependence 
on the contact coupling constant at with respect to S z 
and Si2- A comparison between the two panels shows in 
fact that the characteristic time scale of the oscillations 
in r c is ~ h/a p for S z and S12, and rj ~ h Sij/dp for 
C z , where Sij ~ min{A, g^isB} is the smallest difference 
between eigenvalues of H e . Like in the case of dipolar in- 
teractions discussed above, the leading contributions of 
contact interaction to % are thus quadratic in the hy- 
pcrfinc Hamiltonian for C z , and linear for the other two 
DOF. 

In conclusion, we have shown that the nuclear-induced 
decoherence in a prototypical spin triangle strongly de- 
pends on the qubit encoding. In particular, the decoher- 
ence time of the chirality qubit approaches the ms range, 
i.e. several orders of magnitude larger than the gating 
times predicted for the manipulation by means of elec- 
tric fields. Such robustness results from the decoupling 
of C z from both the total spin orientation and the spin 
texture, preserved by the alignment of the magnetic fields 
in the direction normal to the spin triangle. The eigen- 
states of S12 are characterized instead by different spin 
textures, and thus couple differently to the nuclear spin 
bath. This results in decoherence times comparable to 
those of the total spin projection S z . Larger decoherence 
times for S12 could be obtained in spin clusters where 
the distance between electron spins is much smaller than 
that between electron and nuclear spins. Finally, the 
presence of a strong contact interaction between nuclear 
and electron spins introduces an oscillationg behavior in 
the decoherence factor of all the DOF, with characteris- 
tic time scales that are 10 2 times larger for spin chirality 
than for S z and S12. 

We acknowledge financial support by PRIN of the Ital- 
ian MIUR, by the Swiss NF, and by FP7-ICT project 
"ELFOS". 
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